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1. Introduction 

Bethe-Salpeter (BS) equation for a relativistic bound system was initially formulated in the 
Minkowski space [Q]. It determines the binding energy and the BS amplitude. However, in prac- 
tice, finding the solution in Minkowski space is made difficult due its singular behaviour. The 
singularities are integrable, but the standard approaches for solving integral equation fail. To cir- 
cumvent this problem, the BS equation is usually transformed, by means of the Wick rotation, into 
Euclidean momentum space. 

Some attempts have been recently made to obtain the Minkowski BS amplitudes. The ap- 
proach proposed in [Q] is based on the integral representation of the amplitudes and solutions have 
been obtained for the ladder scalar case [Q, ||] as well as, under some simplifying ansatz, for the 
fermionic one M\. Another approach [f7|] relies on a separable approximation of the kernel which 
leads to analytic solutions. 

In previous works [jlj, ||] we have proposed a new method to find the BS amplitude in Minkowski 
space and applied it to the system of two spinless particles. Like in the papers [Q |5| 0], it is based 
on Nakanishi integral representation of the BS amplitude. The main difference between our 
approach and those followed in [Q, ||, |6|] is the fact that we use the light-front projection. This 
eliminates the singularities related to the BS Minkowski amplitudes. The method is valid for any 
kernel given by the irreducible Feynman graphs. 

We present in this contribution a generalisation to the fermion systems of our preceding works 
§]. A more detailed version can be found in We will see that the direct application to the 
fermionic kernels of the method used in the spinless case, is however married with some numer- 
ical difficulties. Although they could be overcome by a proper treatment of the singularities, in 
this work we propose an alternative method allowing to solve the BS equation for two fermions 
in Minkowski space with the same degree of accuracy than for the scalar case. The numerical 
applications are here limited to the J 71 = + state. 

The main steps in deriving system of equations for the Nakanishi weight functions are ex- 
plained in sect. |2| starting from the original BS equation. Numerical results for the scalar, pseu- 
doscalar and massless vector exchange couplings are presented in sect. ||. Section || contains 
concluding remarks. 



2. System of equations 

We have considered the following fermion Q¥,m) - meson (0,ju) interaction Lagrangians: 
(0 Scalar coupling JS?£ } = gW¥<j>. 
(ii) Pseudoscalar coupling = ig^js^ip. 

(Hi) Massless vector exchange = ^0^. with ILjv = —ig^v/q 2 as vector propagator. 

Each interaction vertex has been regularized with a vertex form factor F(k — k') by the replace- 
ment g — > gF(k — k') and we have chosen F in the form: 

q 2 - A 2 + 1£ 
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Let us first consider the case of the scalar coupling and the corresponding ladder kernel. The 
BS equation for the amplitude <I> reads: 

dV 



i(m + jp + k) 



ikp+kf 



- nf 



IE 



-ig 2 )F 2 (k-k') 



(2k) 



(k-k') 2 -H 2 + i£ 



i(m — Ip + k) 



(\p - k) 2 - m 2 + i£ 



where p = k\ + k 2 , k = {k x - k 2 ) /2, k! = (k[ - k' 2 ) /2. 

In the case of J 71 = + state, the BS amplitude has the following general form: 

&(k,p)=S 1 fa+S 2 fo + S 3 fa+S4fa 



(2.2) 



(2.3) 



where S,- are independent spin structures (4x4 matrices) and (j), are scalar functions of k 2 and p ■ k. 

The choice of 5,- is to some extent arbitrary. To benefit from useful orthogonality properties 
we have taken 



7s, 



where o^ v = 4 {YnYv — 7v7/i)- The antisymmetry of the amplitude (2.3) with respect to the permu- 
tation lo2 implies for the scalar functions: 0i,2,4(&,p) = 0i,2,4( — k,p), ty?,(k,p) = —^(~k,p). 

A decomposition similar to (2.3) was used in ^ to solve the BS equation for a quark-antiquark 
system but the solution was approximated keeping only the first term S\ <pi . 

We substitute (23) in eq. (^2|), multiply it by S; and take traces. As we will see (left panel in 
fig. [I] below), the kernel in the resulted equation, in contrast to the spinless case, is still singular. 
These singularities are integrable. They do not prevent from finding numerical solution, but they 
reduce its precision. This can be avoided by a proper regularization of equation, multiplying both 
sides of it by the factor 

(m 2 -L 2 ) (m 2 -L 2 ) 



n{k,p) 



(2.4) 



[(§ +k) 2 -I? + ie] [{^-k) 2 -L 2 + ie\ 

This factor has the form of a product of two scalar propagators with mass L. It plays the role of form 
factor suppressing the high off-mass shell values of the constituent four-momenta k\ 2 = (f 
and tends to 1 when L — > oo. In this way, we get the following system of equations for the invariant 
functions </>,•: 

r](k,p) 



n(k,p) <Pi(k,p) 



[(§ + k) 2 -m 2 + ie] [(f - k) 2 -m 2 + is] 
d 4 k' ig 2 F 2 (k-k'^ 4 



(2tt) 4 (k-k 



'\2. 



li z + ie 



Y^c ij (k,k / ,p)(j)j(k / ,p), 



(2.5) 



i=i 



Since tj (k,p) ^ 0, the equation thus obtained is strictly equivalent to the original one. We will 
see however that, due to the presence of the r\ factor, the light front projection modifies the resulting 
kernels which become less singular functions. The coefficients c,y are determined by traces and are 
given in [g]. 

Then we represent each of the BS components fy(k,p) by means of the Nakanishi integral 



<pi(k,p) 



-l 



dz' / d/ 



8i(M 



[k 2 +p-kz' + \M 2 -m 2 -Y + ie 



3 ' 



(2.6) 
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z' z 



Figure 1: Left: The kernel matrix elements V\a{j,Z', Y,zf) for z — 0.7, 0.9, 0.95 as a function of z' and fixed 

values of y, /. The discontinuity is finite at a fixed value of z but diverges when z 1. 

Right: Regularized kernel Vf 4 (y,z',Y v.s. z' for fixed values ofz = 0.7, 0.9, 0.95 andL= 1.1m. 



and apply the light-front projection to the set of coupled equations for the corresponding weight 
functions gi(y,z). As mentioned in the Introduction, this projection, which is an essential ingredient 
of our previous works [jlj, [|, 1C], consists in replacing k — > k + in eq. (2.5) and integrating 



over j6 in all the real domain. 

The technical details of the light-front projection are similar to those explained in ref. [§] for 
the case of the spinless particles. We obtain in this way a set of coupled two-dimensional integral 
equations: 

J o d/ J dz' V*{Y,z;Y,z!) gi (Y,z') = Y t J o d/ J Az' Vfjir^Y^gjiY^') (2.7) 

The kernel V 8 and also Vfj for all types of couplings and states are given in |9|]. These kernels 
depend on the parameter L. Closer is L to m, smoother is the kernel and more stable are the 
numerical solutions. However the weight functions gi{y,z) as well as the binding energies provided 



by ( |2.7| ) are independent of L. 

The kernel V g is finite and it vanishes for z = ±1. For a fixed values of y,z and /, V g is a 
continuous function of z' with a discontinuous derivative at z' = z. 



As already mentioned, without Y\ -factor, most of the kernel matrix elements V-f are singular. 



Namely, they are discontinuous at z' = z. In some cases - like e.g. Vu displayed in fig. H (left) - 
the value of the discontinuity, although being finite at fixed value of z, diverges when z — > ±1. 

For T] 7^ 1, i.e., for a finite value of L in ( |2.4[ ), the z'-dependence of the regularized kernels is 
much more smooth and therefore better adapted for obtaining accurate numerical solutions. In fig. 
[T] (right) we plotted the regularized kernel Vf 4 as a function of z' for the same arguments y,z, Y an d 
parameters than in fig. |l| (left), where it was calculated without the Tj (k,p) factor. As one can see, 
the kernel is now a continuous function of z 1 ■ 

We would like to emphasize again that despite the fact that the non-regularized and regular- 
ized kernels are very different from each other (compare e.g. figs. [T], left and right) and that the 
regularized ones strongly depend on the value of L, they provide - up to numerical inaccuracies 
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- the same binding energies and weight functions We construct in this way a family of 

equivalent kernels. 

Table 1: Left: Coupling constant g 2 as a function of binding energy B for the 7 = state with scalar (S), 
pseudoscalar (PS) and massless vector (positronium) exchange kernels. The vertex form factor parameter is 
A = 2 and the parameter of the r\ factor L = 1 . 1 

Right: Coupling constant g 2 BS as a function of binding energy B for the positronium / = state in BS 
equation in the region of stability without vertex form factor (A — > °°), i.e. g < n. They are compared to the 
non relativistic results g^g. 
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3. Numerical results 



The solutions of eq. ( |2.7[ ) have been obtained using the same techniques than in ref [gp. We 
have computed the binding energies, defined as B = 2m — M, and BS amplitudes, for the J = + 
two fermion system interacting with massive scalar (S) and pseudoscalar (PS) exchange kernels and 
for the fermion-antifermion system interacting with massless vector exchange in Feynman gauge. 
In the limit of an infinite vertex form factor parameter A — > °°, the later case would correspond 
to positronium with an arbitrary value of the coupling constant. All the results presented in this 
section are given in the constituent mass units (m = 1) and with L = 1.1. 

For the scalar and pseudoscalar cases, the binding energies obtained with the form factor 
parameter A = 2 are given in the left table ||. We present the results for jj, = 0. 15 and }X = 0.50 boson 
masses. They have been compared to those obtained in a previous calculation in Euclidean space 



[ ]11[ ] using a slightly different form factor. Once taken into account this correction, our scalar results 
are in full agreement (four digits) with [p]]. The pseudoscalar ones show small discrepancies 
(« 0.5%). We have also computed the binding energies by directly solving the fermion BS equation 
the Euclidean space using a method independent of the one used in JTT[]. Our Euclidean results are 
in full agreement with those given in the table |TJ 

The B(g 2 ) dependence for the scalar and pseudoscalar couplings is plotted in figs. ||. Notice 
the different g 2 scales of both dependences. The pseudoscalar binding energies are fast increasing 
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Figure 2: Left: Binding energy for scalar exchange v.s. g 2 for A = 2, L = 1 . 1, jx — 0. 15 and jX = 0.5. 
Right: Binding energy for pseudoscalar exchange v.s. g 2 for A = 2, L = 1 . 1, fi = 0. 15 and ji = 0.5. 



B/m o.: 



Positronium J=0 1 / 




- NR / 




-— BSA=2 /] 




■— BSA=5 / ! 




BS / | 




,1 , , 


i , 1 , , , , 



Figure 3: Binding energy for J=0 positronium state versus g 2 (black solid line) in the stability region g < 
g c = ft- Dashed and dotted-dashed curves correspond to the results for increasing values of the vertex form 
factor parameter A. They are compared to the non relativistic results (red solid line). 



functions of g 2 and thus more sensitive to the accuracy of numerical methods. This sharp behaviour 



K. 



was also exhibit when solving the corresponding light-front equation [12]. 

In the positronium case, we found the existence a critical value of the coupling constant g c 
For g < g c , solutions with finite binding energy exist without form factor (i.e., at A — > «>). 

The ground state positronium binding energies without vertex form factor are given in table |l] 
(right) for values of the coupling below g c , nonrelativistic results gjj R = %%\jBjm are included for 
comparison. One can see that the relativistic effects are repulsive. 

These results are displayed in fig. || (black solid line), and compared to the binding energies 
obtained with two values of the form factor parameter A = 2 (dashed) and A = 5 (dot-dashed). The 
stability region is limited by a vertical dotted line at g = g c = n. Beyond this value the binding 
energy without form factor becomes infinite and we have found B(g — > g c ) ?» 0.10. The inclusion 
of the form form factor has a repulsive effect, i.e. for a fixed value of the coupling constant it 
provides a binding energy of the system which is smaller than in the A — > oo limit (no cut-off). 
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Finally, in figs. || we present some examples of the Nakanishi weigh functions g;(y,z). They 
correspond to a B = 0.1 state with the scalar coupling and the same parameters A = 2, /J. = 0.50 
than in table [[} In the left figure is shown the /-dependence for a fixed value of z and in the right 
figure - the z-dependence for a fixed y. Notice the regular behaviour of these functions as well as 
their well defined parity with respect to z - g\,2A are even an d g3 is odd. 




A=2 (1=0.50 B=0.10 






. . gl -pO.54 




, . g2 T=0.54 




■ g3 -pO.54 




. g4 -pO.54 


" e=10" 6 



Figure 4: Left: Nakanishi weight functions v.s. y for z = 0.6, for scalar exchange for A = 2, L = 1.1, 
ji =0.15 and jj. = 0.5. Right: Nakanishi weight functions v.s. z for 7= 0.54. 

Corresponding BS amplitudes 0; are displayed in figs. |5[ The left panel represents the k de- 
pendence of (f>j for a fixed value of | k |= 0.2. They exhibit a singular behaviour which corresponds 
to the pole of free propagators ko = Ek — y in r.-h.-side of eq. 12. The right panel represents the 
I k I dependence of the amplitudes </>,- for a fixed value ko = 0.04. For this choice of arguments, the 
amplitudes are smooth functions of | k |, though they will be also sing ular for ko > f = 0.05. 




Figure 5: Left: Bethe-Salpeter Minkowski amplitudes, corresponding to fig. 4, v.s. £0 for k = \k\ = 0.2. The 
amplitudes 0i and 02 are indistinguishable. Right: The same as at left v.s. k =| k \ for £o=0-04. 



4. Conclusions 

We have presented a new method to obtain the solutions of the Bethe-Salpeter equation in 
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Minkowski space for the two-fermion system. It is based on a Nakanishi integral representation of 
the amplitude and light-front projection and constitutes a natural extension our previous works for 
the scalar case [0, p. 

A straightforward generalization of this approach, however results into a singular fermionic 
kernel. In order to smooth this singularities, a proper regularization of the kernels has been pro- 
posed. This generates a family of strictly equivalent equations depending on one parameter L. 
Their solution gives the same binding energies and Nakanishi weight functions gi(y,z) . 

The binding energies for the scalar and pseudoscalar exchange kernels and for massless vector 
exchange (positronium) have been found. They coincide with the ones obtained in Euclidean space, 
thus providing a validity test of our method. The solutions for the scalar and positronium states 
without vertex form factor (A — > oo) are found to be stable below the following critical values g c of 
the coupling constant: g c = 2% (scalar) and g c = % (positronium). 

The BS amplitudes in Minkowski are obtained in terms of the computed Nakanishi weight 
functions. They exhibit a singular behaviour due to the poles of the free propagators. 
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